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SUMMARY OF VARIABLE PROPERTY HEAT-TRANSFER EQUATIONS AND 
THEIR APPLICABILITY TO A NUCLEAR ROCKET NOZZLE 

Maynard F. Taylo • 

NASA Lewis Research Center, Cleveland, Ohio 


At act 

'i’h* prediction equal ion resulting from a t horough study of all available 
i iijgle-phase hydrogen heat-transfer data is modified to predict heat-transfer 
coefficients on both the concave and convex surface of curved tubes. The pre- 
diction equations for straight and curved tubes are then applied to the 
Phoebus-2 nuclear rocket nozzle. Predicted exit pressure and temperatures 
are compared with experimental values from recent nuclear tests of the 
Phoebus -2 nozzle, and found to be in very good agreement. 


INTRODUCTION 

The extreme condition, encountered in regeneratively cooled nuclear rocket 
nozzles produce severe heat-transfer problems in the coolant passages. An 
effective method of predicting heat-transfer coefficients in the cooling pas- 
sages is essential to the optimization of any nozzle design. Of particular 
concern is the high heat flux throat region fluxes of 20 Btu per second per 
square inch (32.7 megawatts per square meter) and higher are reached. 

The author of this paper has previously studied all available single-phase 
hydrogen heat transfer for flow through straight tubes and recommended a 
single correlation equation for a wide range of conditions [l]. In the pres- 
ent paprr the calculated heat transfer coefficients are shown to be in very 
close agreement with recent experimental values reported in [2] . More re- 
cently the author has modified the straight tube equation to include the ef- 
fects of curvature from [3]. The modified equation is compared with some 
existing experimental data from single curved tubes [4]. These results and 
recommended applications of the equations have been reported by this author 
[5]. 

Of course, the most severe test of the ability of these equations to calcu- 
late heat transfer in the cooling passages of a nuclear rocket nozzle is to 
compare the predictions with experimental results. In this paper the predic- 
tions are compared with the results obtained from recent nuclear tests of the 
Phoebus -2 nuclear rocket nozzle. 


HEAT-TRANSFER COEFFICIENTS IN STRAIGHT TUBES 

A study of 4622 experimentally determined hydrogen heat-transfer coefficient 
data points [l] resulted in the following prediction equation: 

Nu = 0.023 Re£- 8 Pr£- 4 (l) 


C 1 --- (0.57 


1.59\ 
x/D~ ) 


where 


2 


for the range of inlet pressure and inlet temperature shown in figure 1. The 
use of the temperature -entropy diagram is a convenient method of showing the 
location of measured inlet pressure and temperature in relation to the satura- 
tion lines and critical pressure p c and critical temperature T c . 

The hydrogen data are separated into regions 1 and 2 as shown in figure 1. 
These regions are the result of the study of the effect of inlet temperature 
, inlet pressure p^, and the transposed critical temperature (the tempera- 
ture at which the specific heat at constant pressure reaches a maximum) T* . 

In region 1, 87 percent of the 3674 heat-transfer coefficients predicted by 
equation (l) deviates less than ±25 percent from the measured values. Re- 
gion 2 is defined by 45° R < T* < T* and Pc < Pi < 530 psia and is often 
referred to as the near-critical region. In region 2, 40 percent of the 948 
predicted heat-transfer coefficients deviates less than ±25 percent from the 
measured values. At present, there is considerable doubt about the transport 
properties in this region. 

Equation (l) is used to predict heat-transfer coefficients which can be com- 
pared to the most recent experimental data for the range of conditions en- 
countered in nuclear rocket nozzles. The ratio of these predicted coeffi- 
cients to the measured values is shown as a function of the ratio of wall to 
fluid bulk temperatures in figure 2. Figure 2 shows that 94 percent of the 
heat-transfer coefficients predicted using equation (l) deviated less than 
±20 percent from the measured values. 


HEAT-TRANSFER COEFFICIENTS IN SYMMETRICALLY HEATED CURVED TUBES 


Heat -transfer measurements have been reported for both the concave or "swept" 
side and the convex or "unswept" side of curved circular tubes [4], Equa- 
tion (l) was first tried without any change from the form that correlated 
straight tube data to predict heat-transfer coefficients for both the concave 
and convex sides of the symmetrically heated test sections. The rat -0 of the 
heat-transfer coefficient predicted by equation (l) to the measured coeffi- 
cient is shown as a function of temperature ratic in figure 3. As might be 
expected, the predicted heat-transfer coefficient is lower than the measured 
value on the concave side and higher than the measured value on the convex 
side. The calculated heat-transfer coefficient is as small as one-half the 
experimental value on the concave side and as large as almost twice the ex- 
perimental value on the convex side. 

Equation (l) is modified [5] with the Ito correction factor for curvature 
[3] to give the equation 

0.8 0.4/% Y Cl r rfl 0-05 

Nu b = 0.023 Re,, Pr b j |Re b ^)J (2) 


for the concave side. 

The convex side of a curved tube is rarely given any consideration, but in 
reference [5] the reciprocal of the Ito factor was used to modify equa- 
tion (l) to give 
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Nu b * 0.023 Reg- 8 Prg- 4 (^) 1 (3) 

for the convex Bide. 

The results of using equations (2) and (3) to predict the heat-transfer coef- 
ficients on the concave and convex sides of a tube are shown in figures 4(a) 
and (b), respectively. The experimental data in figure 4 is the same data 
used in figure 3. Approximately 90 percent of the predicted heat-transfer 
coefficients deviated less than ±20 percent from the experimental values. 

It appears that the Ito correction does improve the correlation and that 
equations (2) and (3) will predict heat-transfer coefficients on the concave 
and convex sides of curved tubes, respectively. 


ASYMMETRICALLY HEATED CURVED TUBES 


Both the straight and curved tube data presented up to this point are for 
symmetrically heated circular tubes. In the rocket nozzle the cooling pas- 
sages are noncircular and they are heated from only one side. Some single 
tube data are available for asymmetrically heated circular and noncircular 
passages curved to simulate the throat and high flux region of the Phoebus -2 
nozzle [2]. Heat-transfer coefficients predicted by equation (2) were com- 
pared to these experimental values for conditions as near as possible to those 
for the actual nozzle. 

Figure 5 shows the ratio of the heat -transfer coefficients calculated by 
equation (2) to the experimental coefficients as a function of the dimen- 
sionless distance from the entrances. Figure 5 shows the ratio of coeffi- 
cients for the Phoebus -2 contour. The agreement between predicted and mea- 
sured heat-transfer coefficients is good for single, asymmetrical, curved 
tubes that simulate the Phoebus -2 geometry. 

An important use of the heat-transfer coefficient is in the prediction of 
wall temperatures. Equation (2) is used to predict the inside wall temper- 
atures which are compared with the measured wall temperatures [2] in fig- 
ure 6 for the simulated Phoebus -2 coolant tube. The agreement is good and 
the prediction is conservative which is desirable. 


APPLICATION OF PREDICTION EQUATIONS TO ROCKET NOZZLES 


Regeneratively cooled rocket nozzles are made up of noncircular passages 
formed to give the desired area ratios, and are essentially combinations of 
straight and curved tubes. The success demonstrated in predicting heat- 
transfer coefficients in straight tubes using equation (l) and curved tubes 
using equations (2) and (3) encourages the use of these equations in pre- 
dicting heat transfer coefficients in the coolant passages of a rocket nozzle. 

To predict the heat-transfer coefficients in the cooling passage of a nozzle, 
the following equations are recommended [5]. 
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tot ranee of Coolant Passage 


Straight entrance x/P > 1 

Nu b » 0.023 Re b ‘ 0 Pr b * 4 

J0° angle bend and orifice entrance x/P > 5 



Nu b « 0.023 Re b #[ 'Pr b * 4 



( 1 ) 


(4) 


where = 5 for e 90° angle bend and 11 for an orifice entrance. The 
terra (l + F-^D/x)) has been previously reported [6]. 


Throat section (concave curvature) 

/T \ - c i r g "| 0,05 

Nu b - 0.023 Reg- 8 Prg >4 (r^ g) J 

Exit section (convex curvature) 

Nu b - 0.023 Reg- a Pr°' 4 (^j 1 Jr^ 

Any straight sections 



0.8„ 0.4 


0.023 ReJ Pr b 



(2) 


(3) 


( 1 ) 


At present it is not known how far the effects of curvature on the heat- 
transfer coefficient will extend downstream of the point of tangency where 
a curved tube becomes straight. It is reasonable to assume that the effect 
of curvature would diminish with x/D rather than change abruptly at the 
point of tangency. The curvature effect appears to be present at the last 
instrumented station shown in figure 6(b) which is 2 diameters downstream 
of the tangency point. In the calculations presented in this paper the ef- 
fects of curvature was assumed to begin or end at the point of tangency with 
the straight sections. 


COMPARISON OF PREDICTED RESULTS WITH PHOEBUS-2 NUCLEAR TESTS 


A severe test of the prediction equations proposed in this investigation 
is to use them to predict the heat-transfer in the coolant passages of the 
Phoebus-2 nuclear rocket nozzle. The Phoebus-2 nozzle was run at the highest 
power level ever achieved by a nuclear rocket engine, 4080 megawatts. The 
conditions for a wide range of tests are shown in table I. The equations 
presented herein were used to revise the coolant side heat-transfer calcula- 
tions in an existing digital computer program [7] for calculating heat trans- 
fer and fluid flow for convectively cooled rocket nozzles. 

The calculation of heat transfer from the hot propellant is also necessary 
and the Nusselt equation: 



5 


m no 0.3 fc\ 

NUf « ^g Re f ^ r f (-) 

is used in the computer program. Common practice [7] is to use a gas '•r-ef- 
ficient, Cg, that is a function of r.ozzle area ratio. In this investigation 
both a constant C g of 0.026 and an empirically determined variable C g 
provided by the Aerojet General Corp. are used and the results compared. The 
variable C g and the area ratio for the 34 stations in the nozzle is shown 
in table II. There is no provision in the computer program for the small 
amount of film cooling on the hot-gas side of the chamber wall. The maximum 
effect can be estimated by assuming no heat addition at the last 4 stations 
(about 12.5 in. or 32 cm). 

X 5 th friction and momentum pressure drops in the coolant passages were calcu- 
lated in the computer program. The equation used to calculate these pressure 

drops in the program were not changed. The pressure drop is affected by the 

coolant temperature, thus making temperature rise and pressure drop inter- 
related. 

Total temperature and static pressure of the coolant at each axia] station 
for test EP IV-C are calculated and shown in figure 7. The only measured 
values for comparison art the coolant passage exit temperature and pres- 
sure, which are also shown in figure 7. The solid line represents the calcu- 
lated values with C g varying as shown in table III. The dashed line re- 
presents the calculated values using a constant Cg of 0.026. The agreement 

between both calculated values and the measured values is very good. The 

measured and calculated exit pressures and temperatures and percent deviation 
for several other EP-IV tests are shown in table III. For the case of no 
heat addition in the lest four stations (corresponding to film cooling of the 
chamber wall) both the exit temperature and pressure fell within the range of 
accuracy of the measured values and are not shown in table III. 

In figure Q the local hot-gas side wall temperature and the local coolant 
side heat flux is shown as a function of axial position. The solid line 
represents the values calculated using the variable Cg from table III and 
the dashed line represents the values resulting from tne use of a constant 
Cg of 0.026. The greatest differences resulting from the calculations 
using the different C g 's are in the chamber. The use of temperature rise 
and pressure drop measurements indicates little of what happens locally in 
a rc-ket nozzle. Wall temperature and heat flux did not vary as much in 
the throat region as they did in the nozzle chamber. This is expected since 
the Cg in table II is two to three times the constant value of 0.026. 

The good agreement between the two sets of calculations and the measured 
exit temperature and pressure indicates that with accurate predictions on 
the coolant side the hot-gas side equations are of secondary importance in 
predicting exit conditions. However, as one might expect, local wall tem- 
peratures and heat flux are very much affected by local gas -side heat- 
transfer coefficients. 

The rather abrupt changes in wall temperature is due to the application of 
the curvature correction which increased the heat-transfer coefficient on 
the concave or swept surface (decreasing the wall temperature) and de- 
creased the heat-transfer coefficient on the convex or unswept surface (in- 

creasing the wall temperature). In the actual nozzle wall these changes 
might not be so abrupt because of axial heat conduction in the coolant pas- 
sage wall . Axial heat conduction was not accounted for in the computer 

program. 
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CONCLUDING REMARKS 


The computer program for calculating heat transfer and fluid flow In reger- 
<ratively cooled rocket nozzles [7] was revised to include the heaL -transfer 
prediction equations recommended in this invest lgation. 'Hie revised p gran, 
was thi n used to predict exit temperatures and pressures, nozzle wall tem- 
peratures, and heat flux to the coolant. Only the exit temperature and 
pressures from the Phoebus-2 nuclear tests were available for comparison with 
predicted values. Two setr of calculations were made for each nuclear test: 
one in which the gas constant C^. was constant at 0.026, and another test 
which used recommended C R *s that, varied from 0.018 to 0.080. Both sets of cal- 
culations yielded exit temperatures and pressures that are in good agreement 
ith measured values. 'This indicates that with accurate predictions on the 
coolant side, the less well developed heat -transfer correlations on the hot- 
gas side are of secondary importance in predicting total temperature rise and 
static pressure drop. The hot-gas side heat-transfer does, however, change 
the wall tempt rature distribution. This indicates that measured local wall 
temperatures are needed to verify any hot-gaj side equations. The good 
agreement between measured and calculated exit temperature and pressure does 
not necessarily mean that all the heat-transfer and fluid flow problems of 
regeneratively cooled rocket engines are now solved. Instead, it indicates 
a need for work on the hot-gas side and a need for local wall temperature 
and coolant pressure measurements. Such measurements for any rocket engine 
test, especially a nuclear engine, presents a considerable challenge to 
present technology. 


NOMENCLATURE 

C coefficient in Nusselt equation 

Op specific heat of gas at constant 
pressure 

D inside diameter of test section 

G mass flow rate per unit cross- 
sectional area 

h local heat-transfer coefficient 

k thermal conductivity of gas 

Nu Nusselt number, hD/k 

Pr Prandtl number, Cpu/k 

p absolute static pressure 

q rate of heat transfer to gas 

per unit area 

R radius of curvature 

Re^ bulk Reynolds number, GD/u^ 

Rep modified film Reynolds number, 
Pf v b D /^f 

r inside radius of passage 

T temperature 


T* transposed critical temperature 
(temperature at which specific 
heat of fluid at constant 
pressure reaches a maximum) 

v velocity 

x distance from entrance of test 
section 

u absolute viscosity of gas 

P density of gas 

Subscripts : 

b bulk (when applied to properties, 

indicates evaluation at bulk 
temperature T^ ) 

c critical 

f film (when applied to properties 

indicates evaluation at film 
temperature Tf) 

g gas -side 

i inlet 

w wall 
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TABLE I - OPERATING CONDITIONS FOB PHOEBUS 2A EP-IV NUCLEAR ROCKET IT STS * 


Ui U. S Customary Units 


Test 

Flow rate, 
lbm /sec 

Chandler condition* 

No//|e inlet conditions 

Thermal power, 
MW 

Hot pas 

Coolant 

Pressure, 

psla 

Teni|ie rature, 

"r 

Pressure. 

psla 

Temperature. 

°R 


10H 6 

97, 9 

124 


197 

40 8 




125 6 

214 


362 

43.0 



192 2 

176 6 

351 


650 

47. 2 



244 1 

219 8 

435 

2940 

820 

50. 9 



245. 2 

219 5 

506 

3890 

975 

52 1 


6 

262 3 

234 9 

555 

4060 

1080 

53 9 

4080 


(l)l S I. Units 


Test 

Flow rate, 
kp sec 

Chandler conditions 

Nozzle inlet conditions 

Thermal power, 
MW 

Hot pas 

Coolant 

Pressure, 
MN m 2 

Temperature. 

K 

Pressure, 
MN / m 2 

Temperature . 
K 

n 

49 3 

44 4 

0. 855 

689 

1. 36 

22. 6 

490 

n 

64 4 

57 0 

1 48 

1178 

2. 50 

23 9 

1090 


87. 2 

80 l 

2. 42 

1711 

4. 48 

26 2 

2190 

n 

110 7 

99. 7 

3. 00 

1633 

5. 65 

28 3 


Rfl 

111. 2 

99. 6 

3 49 

2161 

6. 72 

28. 9 


D 

119.0 

106 6 

3 83 

2256 

7.45 

29. 9 

4080 


‘‘Unpublished data received from Donald L Hanson, Los Alamos Scientific Laboratory. 
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TABLE II - AH FA HATH) AND GAS COEFFICIENT 


FOR EACH CALCULATION STATION* 1 


Region 

Station 

Area 

ratio 

Gas 

coefficient, 

c « 

Region 



C.as 

coellicient, 

C 

g 

Divergent 

1 

7 00 

wmm 

Convergent 

21 

2 

34 

0. 027 



6 44 

:flfCfr 


22 

3. 

17 

.EJ. 



5. 90 

.033 


23 

4. 

13 

EX 



5 39 

.033 


24 

5. 

21 

.029 



4 90 

034 


25 

6.42 

032 


6 

4 43 




26 

7. 

76 

. 03 !» 


7 

3. 98 




27 

9 22 

. 038 


8 

3 56 






— 



9 

3. 17 



Knuckle 

28 

10. 81 

0 042 


10 

2 79 




29 

11. 85 

. 045 


11 

2 44 




30 


22 

. 048 


12 

2 12 

. 032 

Chamber 

31 

12. 22 



13 

1.81 

. 030 


32 



. 064 


14 

1 53 

028 


33 



. 075 


15 

1. 28 

.024 


34 



. 080 

Throat 

16 

1.04 








17 

1. 00 

ESI 







16 

1 07 

D 







19 

1. 27 

Cl 







20 

1. 63 

.029 







‘‘Unpublished data received from .lames O. Sane, Aerojet-General Corp. 
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E-5221 



0 10 20 30 40 50 

Inlet entropy, Jf(gK°KI 


Figure 1. - Range of hydrogen inlet temperature and inlet pressure for which equa 
tion ill has been experimentally checked. 


Figure 2. * Variation of the ratio of calculated 
to measured heat transfer coefficients witt 
temperature ratio. Straight tube data [2]. 





Calculated h/Measured h 


<ai Concave. 



<b) Convex surface, equation (1), 

Figure 3. - Variation of the ratio of calculated 
to measured heat transfer coefficients with 
temperature ratio. Data from symmetrically 
heated curved tubes [4], 



(b) Convex surface, equation (3). 

Figure 4. - Variation of the ratio of calculated 
to measured heat transfer coefficients with 
temperature ratio. Data from symmetrically 
heated curved tubes [4]. 
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Assymetrical Phoebus*2 simulated tube, h calculated 
using equation 12) 


Figure 5. - Variation of the ratio of calculated to measured heat 
transfer coefficients with axial distance from entrance. 

Data from asymmetrically heated curved tube [2], 



Figure 6. - Variation of inside wall temperature with axial 
distance from entrance for the Phoebus-2 contour [ZJ. 
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Figure 7. - Variation of coolant total tempera- 
ture and static pressure with axial location. 
Thermal power 4060 megawatts. 



